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Abstract 



We point out that the 1/A^ expansion, which is widely invoked to infer 

Q^' properties of the 2D 0{N) models, is nonuniform in the temperature, 

D . i.e. with decreasing temperature the 1/A^ expansion truncated at a fixed 

order deviates more and more from the true answer. This fact precludes 

the use of the expansion to deduce low temperature properties such as 

^ . asymptotic scaling for those models. By contrast, in the ID 0{N) 

H ' chains, there are no signs of such a nonuniformity. 
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1 Introduction 

The 1/A^ expansion [|l| is a very popular tool to analyze 2D 0{N) models, 
especially in the critical dimension 2. But already the first rigorous paper 
by A.Kupiainen p[, which established the fact that the expansion is a valid 
asymptotic expansion in the thermodynamic limit at fixed (3 = jS/N, contained 
hints that the expansion is really only good at high enough temperature: the 
main result stated that the expansion is valid provided (3 is less than the critical 
value of the 'spherical model', i.e. the N —>■ oo limit of the model. Likewise the 
foUowup paper by Frohlich et al [§] proved Borel summability of the expansion 
only for small /?. Miiller, Raddatz and Riihl |^ computed the 1/A^ corrections 
to the mass gap and the susceptibility for large /3; their results are nonuniform 
in (3 in the sense that the corrections grow beyond all bounds as /3 —>■ cxd. 
This fact, which is apparent from the equations in was, however, not noted 
explicitely by the authors. Finally, U.Wolff computed the correlation length 
of the 0{N) model for A^ = 2, 3, 4, 8. Comparing his results with the spherical 
model makes it quite obvious that the approach to the A^ ^ oo limit slows 
down with increasing (3 (see Tab.l). 

In spite of all those facts, many authors (e.g. U], H, |^) continue to cite the 
1/N expansion and the properties of the spherical model as evidence for the 
supposed asymptotic scaling and absence of a critical point of the 2D 0{N) 
model at finite A^. In this note we report a few more calculations that show 
explicitly that the 0{1/N) corrections for a number of quantities grow linearly 
with f3. All those quantities are of long range character, such as correlation 
lengths, susceptibilities or two-point functions at large distances. We did not 
see such effects in short range quantities such as the energy (= two point 
function at one lattice distance); this does not mean, however, that the 1/A^ 
expansion for short range quantities is uniform in (3. Uniformity would require 
control of the deviation of the true result from its truncated 1/A^ expansion, 
something that is not analyzed here. We should, however, point out that our 
recent results on the problems of perturbation theory (PT) ^ in non-Abelian 
models suggest that in fact there is nonuniformity even in those short- distance 
quantitities. 

The situation is different in ID: no indication of any non-uniformity is 
found. The underlying reason is that in ID the necessary computations can 
be reduced to the analysis of simple 1 dimensional integrals, so that no infrared 
effects, which in our mind are responsible for the nonuniformity in 2D, can 
arise. 

2 Notation and Definitions 

Consider the lattice 7L and let A be some given finite subset of it. The non- 
linear 0{N) o"-model with standard nearest neighbor interaction (s.n.n.i.) at 



inverse temperature (3 = (3N is defined by the generating function 

Zk{.g) = / Ylexpi^pSi ■ Sj)'[[(exp{gi ■ Si)(5(s- - l)o?sJ (1) 

(ij) i ^ ^ 

or equivalently tlie Gibbs measure 

-i- nexp(/?.. ■ s,)X{{5{sj - l)ds) (2) 

Here s is a A^-vector, tlie integration is over tlie spins inside A witli some 
boundary conditions. 

It is well known (see for instance 0) how one can develop an asymptotic 
expansion in powers of 1/N by introducing a 'Langrange multiplier field' a 
dual to the ^-function constraint and a standard saddle point expansion. In 
terms of the Lagrange multiplier field one obtains a dual representation of the 
model; the generating function is now given by 

Z{g)=jX{da,eM-NS)ex^[^-^{g,{-^ + m'-^r^g)^ (3) 

where the dual action is given by 

S = -tr ln(-A + m^ - ^) - z/3 V 4^ (4) 

The mass parameter m appearing in this equation is fixed by the gap equation 

/3 = (-A + m')-^ (5) 

Kupiainen showed, using the technique of cluster expansions, that the re- 
sulting asymptotic expansion is uniform in the volume A and can be used 
to obtain a valid asymptotic 1/A^ expansion in the infinite volume, provided 
(3 < Psph, where Psph is the critical point of the spherical model. It is well 
known how the 1/A^ expansion can be cast in the form of Feynman graphs 
involving two types of propagators B and C given by 



B-,' = -iCi) (7) 



Q, = (-A + m%' (6) 

1 

Those Feynman graphs can be easily evaluated numerically on finite lat- 
tices. In this paper we are considering the 1/A^ corrections to the correlation 
length ^, the susceptibility x and the following ratio of two-point functions 

G(W2) 



where .^oo is the correlation length of the spherical model (up to a correction 
that is negligible for (3 > .5, ^oo = 1/itl where m is the solution of the gap 
equation) and 

G{x) = (s(0) ■ s{x)) (9) 

The quantities ^, x^ ^ ^^ w^ll as the energy E are expanded as follows: 

e = Ui + ^ei + o(^)) (10) 

X = Xoo(l + ^Xi + 0(^)); (11) 

§§||^A.(l + >. + 0(i,)) (12) 

E ^ (s(0) ■ s{l)) = E^l + lei + O(^) (13) 

In the following we will study the behavior of the coefficients .^i, Xi ^^^ '^i for 
large /?. 



3 What is the problem? 

Let us consider a certain physical quantity A{N, 13). According to the preceding 
discussion it will have an asymptotic expansion of the form 

fc=0 

The meaning of the symbol ~ is that the series on the right hand side, if 
truncated at the term of order N~^., will approximate the left hand side by up 
to a remainder term which is o{N~k), and this is true for any k > 0. That 
means 

k 



r,iNj)^AiNj)-j:^ (15) 

fc=0 

satisfies 

lim \rk{N,(3)\N'' = (16) 

N^oo 

Uniform asymtoticity means the stronger property 

\rk{Nj)\N''<Ck{N) (17) 

where Ck{N) is independent of j3 and satisfies 

lim Cfc = (18) 

Af^oo 



It is easy to see that uniform asymtopticity does not hold if any one of the 
expansion coefficients ak{P) is unbounded for large /3, but it should be stressed 
that this is not a necessary condition. There are simple examples of functions 
where every one of the coefficients Cfc(/3) is bounded and yet the expansion is 
not uniformly asymptotic. Consider for instance the function 

f{NP)=PeM-N/(3) (19) 

Its asymptotic expansion in powers of 1/A^ has all coefficients ak{(3) = and 
the remainder term is always equal to /, which is unbounded. 

4 One Dimension 

In ID the situation is particularly simple because the thermodynamic limit of 
the expectation values of local quantities can be obtained by computing in a 
finite box with free boundary conditions. For nearest neighbor quantities like 
the energy E one only has to evaluate a ID integral; long range quantities 
such as the correlation length ^ or the susceptibility x can be expressed in 
terms of E. It turns out that for this reason no nonuniformity in any of those 
quantities is found. 

The energy of a ID 0{N) chain can be obtained as follows: 

E^{s{0)-s{l)) = ^\nZ (20) 

where Z is the partition function for one link: 

' dte^'il-e)"^ (21) 

-1 

Instead of using the general scheme outlined in the previous section, it is 
convenient to obtain the 1/A^ expansion directly from this expression. The 
two point function in the thermodynamic limit is simply given by 

(s(0) ■ s{n)) = E- (22) 

and hence the correlation length is 

C = -\nE (23) 

and the susceptibility 

Introducing the rescaled inverse temperature 



''=N^. P^) 



the partition function takes the form 



' f/te(^-3)/(*) (26) 

-1 

with 

f{t)=^t+^\n{l-t') (27) 

It is now straightforward to produce a saddle point expansion in powers of 
1/(A^ — 3): the saddle point value to is given by 



and the logarithm of the partition function is expanded as 

1„Z = (^-3)/,,,,, + 1iu^^^-^^;L_ + o(_1_) (29) 

Since to is a function that increases monotonically from to 1 as /3 increases 
from to oo, we may use to instead of (3 as an independent variable. Using 



d 1-tiY d 



df3 l + tl dtc 
the expansion for the energy becomes 



(30) 



to il-tl){3 + 2tl) 1 

^-'°"iv^ TT^^ + ^^(iv33F^ ^^^ 

It is already visible that the ration of the 1/(A^ — 3) correction to the leading 
term remains bounded as (3 —>■ oo, i.e. to -^ 1; in fact it is 0(1/ f3). But we 
want to expand in powers of 1/A^, not 1/(A^ — 3), hence we have to re-expand 
the leading term, using 

Again it is very easy to see that we just obtain another correction of order l/j3 
to the 1/A^ contribution ei (see eq.(13)). It is also easy to see that the result 
is in agreement with the one found by Hasenfratz |Q for the infinite chain. 

Using this information and the formulae above it is now trivial to conclude 
that 



261^2 

^^ = -{1 + E){1-E) ^'"^ 

Oi = — ei Ini? (35) 



are all 0(1) for /3 ^ oo as claimed. 



5 Two Dimensions 

We computed the Feynman graphs for the 1/A^ correction to the two point 
function numerically on lattices satisfying L/S,oo = 7. It is observed that this 
is sufficient to reach the thermodynamic limit with a precision of 1 to 2%. 
The expansion for the two point function has the form 

G{^) = Goo(^) + ^ E ^^oo(^ - jMj - 0Goo(0 (36) 

where Gooii) = C{i)/P and the 'self-energy' S is given by 



^ j,l,k€A ^ 

(37) 
and analogously in momentum space 

G{v) = Goo{p)'Hp) (38) 

where S is given by 



Mp) = - E ^Goo{k)B{P -k) + Ib{0) E ^,Goo{k)'G^{k + q)B{q) (39) 

k k,q 

the momenta p,k,q range over 27ri/L, i = 1,2,...L. From the momentum 
space two-point functions G(0) and G{2ti/ L) we obtain the susceptibility and 
also the (effective) correlation length using the formula 



e 



2sinf^G(^ 



^^°) 1 (40) 



which is expanded in powers of 1/A^ using the expansion for the two-point 
functions. 



The results are given in Table 2. They show clearly that the 1/A^ corrections 
for all three quantities grow linearly with /5, respectively Inm, in agreement 
with the findings of ^. 

6 Discussion 

We have found by explicit calculation that the 1/A^ expansion is not uniformly 
asymptotic in /3 for a variety of long range quantities in the 2D 0{N) model. 
The nonuniformity appears in the most blatant form: the correction grows 
linearly with /3. For short distance quantities like the energy we could not find 
this kind of phenomenon up to 0{1/N^). This does not mean, however, that 
the 1/A^ expansion in uniform for those quantities; to show uniformity one 
would have to control the remainder terms. 

This situation is in striking contrast with the one for the ID chains: since all 
compututations reduce to one-dimensional integrals, there is no nonuniformity 
occurring. 

Some authors (for instance 0) try to improve the quality of the 1/A^ ex- 
pansion by introducing a certain 'mass renormalization'. Since we are working 
on a lattice, there is first of all no need to do so since there are no divergencies 
to be absorbed. Furthermore, the procedure corresponds to a a resummation 
of certain terms in the 1/A^ expansion and therefore is no longer a system- 
atic asymptotic expansion in any parameter. While apparently improving the 
agreement with the Monte Carlo data in the (3 region studied, the mathemat- 
ical status of the procedure remains unclear and it can be expected to reveal 
the same problems as the systematic expansion at larger jS. 

In PI we found that PT for the energy yields results that depend on the 
boundary conditions to order l//3^; this dependence survives the limit N — »• oo. 
In our mind this is a strong indication that the 1/A^ expansion is nonuniform 
in (3 even for short distance quantities, since the spherical limit of the model 
does not show this phenomenon. 

But the main conclusion is that that the nonuniformity for the mass gap 
and similar quantities implies that it also holds for the Callan-Symanzik f3- 
function. This means that the 1/A^ expansion cannot be used to draw any 
conclusions about the phase structure of the model for finite N, in particular 
the fact that the spherical model does not have a soft phase in 2D is in no way 
in confiict with the arguments advanced by us |^, |T^ that strongly suggest the 
existence of such a phase for any finite A^. 



References 



[1] M.Campostrini and P.Rossi, Riv.Nuovo Cini. 16 (1993) 1 
[2] A.Kupiainen, Commun.Math.Phys. 73 (1980) 273 
[3] P.Hasenfratz, Phys.Lett. 141B (1984) 385 

[4] J.Frohlich, A.Mardin and V.Rivasseau, Commun.Math.Phys. 86 

(1982) 87 

[5] V.F.Mliller, T.Raddatz and W.Riihl, Nucl.Phys. B251[FS13] 
(1985) 212 

[6] H.Flyvbjerg and S.Varsted, Nucl.Phys. B344 (1990) 646; 
H.Flyvbjerg and F.Larsen, Phys.Lett. B266 (1991) 92; H.Flyvbjerg 
and F.Larsen, Phys.Lett. B266 (1991) 99; 

[7] P.Biscari, M.Campostrini and P.Rossi, Phys.Lett. B242] (1990) 
225 

[8] A.Patrascioiu and E.Seiler, Superinstantons and the Reliability of 
Perturbation Theory in Non-Abelian Mode/s preprint MPI-Ph/93-87 
and AZPH-TH/93-33 

[9] A. Patrascioiu, AZPH-TH/91-49 

[10] A.Patrascioiu and E.Seiler, iVuci.Piiys.fProc.Suppi. B30 (1993) 184 

[11] R.Gupta and C.F.Baillie, Phys.Rev. B45 (1992) 2883 

[12] J.Apostolakis, C.F.Baillie and G.F.Fox, Phys.Rev. D43 (1991) 
2687 

[13] R.G.Edwards, S.J.Ferreira, J. Goodman and A.D.Sokal, Nucl.Phys. 
B380 (1992) 621 

[14] U.Wolff, Phys.Lett. B248 (1990) 335 



Table 1: Comparison of Monte Carlo results with the spherical limit; Monte 



Carlo data taken from 11(0(2)), [g (0(3)), 
0(8)). The numbers represent ^/Coo- 



13 (0(4)) andH (0(4) and 



p 


.500 


.5667 


.575 


.600 


.650 


.725 


N = 8 


1.307 




1.494 




1.711 


1.983 


N = A 


1.866 




2.857 


3.301 


4.425 




N = 3 


2.644 


5.479 


8.348 


15.629 






N = 2 


9.691 


00 


00 


00 


00 


CX3 



Table 2: The 1/N corrections to various long range quantities. For the defi- 
nition of the quantities given see eq. (10), (11), (12). 



(3 


.49349 


.60600 


.67105 


.71704 


.78173 


.82758 


L 


28 


56 


84 


112 


168 


224 


eoo 


4 


8 


12 


16 


24 


32 


ei 


1.764 


2.927 


3.640 


4.152 


4.879 


5.397 


Xi 


2.860 


4.932 


6.229 


7.171 


8.518 


9.487 


ai 


-.486 


-.853 


-1.058 


-1.205 


-1.416 


-1.567 



